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Abstract 

Dynamical parton densities, generated radiatively from valence-like inputs at 
some low resolution scale, are confronted with recent small-x data on deep inelastic 
and other hard scattering processes. It is shown that within theoretical uncertain- 
ties our previous (1994) dynamical/radiative parton distributions are compatible 
with most recent data and still applicable within the restricted accuracy margins 
of the presently available next-to-leading order calculations. Due to recent high 
precision measurements we also present an updated, more accurate, version of our 
(valence-like) dynamical input distributions. Furthermore, our perturbatively sta- 
ble parameter-free dynamical predictions are extended to the extremely small-x 
region, 10~ 8 < x < 10~ 5 , relevant to questions concerning ultra-high-energy cosmic 
ray and neutrino astronomy. 



1 Introduction 



The guiding physical idea underlying the dynamical (radiative) parton model is that 
the steep behavior of the momentum distributions xf(x,Q 2 ) (/ = q, q, g) at x < 10~ 2 , 
Q 2 > 1 GeV 2 is a purely perturbative phenomenon. In fact, in [| non-steep (valence- 
like) initial distributions xf(x,fi 2 ) at some low scale \i « 0.6 GeV were suggested in 
order to predict, purely dynamically, the rise of xq and xg in the above range of x and 
Q 2 . This prediction was subsequently confirmed by the measurements of F%{x, Q 2 ) and 
xg(x,Q 2 ) at HERA |3], ||. As was stated in [^, ||, the available pre-HERA data at 
x > 10~ 2 utilized to fix the valence-like input distributions still allowed for a slight, 
typically 0(10%), variation (increase) of \i which did not affect F%{x > 10~ 2 , Q 2 ) but 
resulted in an about 10% (20%) uncertainty of the radiative predictions at x = 10~ 3 (10~ 4 ). 
Hence the O(10%) discrepancies between the distributions of and recent precision 
measurements at HERA || [F], |8| are obviously neither unexpected, nor do they invalidate 
the notion of a radiative, i.e. dynamical, origin for the steep rise at x < 10~ 2 and Q 2 > 
1 GeV 2 . Indeed, a fine-tuning of \i and/or /(x,/i 2 ) was always understood |2|, [3| to be 
necessary in due course. 

Given the accuracy of recent HERA small-x data, as well as new large-x constraints, 
mainly on the flavor decomposition of the quark sea, it is now appropriate to perform an 
update of our previous dynamical parton distributions and to follow the effects of these 
fine-tunings on the predictions for F 2 (x,Q 2 ), g(x,Q 2 ) and other relevant deep inelastic 
observables. 

Section 2 will be devoted mainly to a discussion concerning the small-x implications 
of the recent HI and ZEUS high precision data |], [7], [| on F^x.Q 2 ) and g(x,Q 2 ) and 
how they (slightly) modify our previous valence-like LO/NLO input gluon and sea 
distributions. Further data on F£ (x, Q 2 ) and fixed-target results on Ff , F 2 n / and d v / u v , 
as well as data on pp(d) — > fi + [i~X and pp — > W ± X asymmetries relevant for fixing and 
testing our g, u v , d v , d — u and u/d densities are presented and compared with the present 
update and with our previous GRV(94) results in Section 3. Furthermore, we also extend 
here our perturbatively stable and parameter-free dynamical small-x predictions to the 
ultra small-x region (10 -8 < x < 10~ 5 ), relevant to questions concerning ultra-high-energy 
cosmic ray and neutrino astronomy. Finally, in Section 4 we present a brief summary and 
a general discussion concerning the status of the dynamical (radiative) parton model in 
the light of present and future data and its application in theoretical NLO analyses. 
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2 Consequences of recent data on F 2 (x, Q 2 ) and parton 
densities 



As stated in the Introduction we intend to study the implications of modifications of our 
original GRV(94) input ||. We shall proceed in a stepwise manner and begin with a 
mere modification of the input scale /z 2 keeping everything else unchanged. The results 
shown in Figs. 1 and 2 illustrate the point that the previously noted B [| uncertainty in 
/i 2 can accommodate the main discrepancies between || and the recent data || [7|, |j on 
F 2 (x, Q 2 ) and the experimentally extracted xg(x, Q 2 ). It should be noted that our present 
new (modified) analysis yields results which lie in between the two curves shown in Fig. 1 
as will be discussed and illustrated in Fig. 3. Our new results (solid curve) for xg(x, Q 2 ) 
at Q 2 = 20 GeV 2 are shown as well in Fig. 2 where they practically coincide with the 
one of GRV(94) with /^nlo increased by 10% (dotted curve). The fact that the somewhat 
too steep GRV(94) predictions for F 2 (x,Q 2 ) and xg(x,Q 2 ) at, say, x = 10~ 4 can both 
be corrected with the new slightly modified valence-like input constitutes a nontrivial 
confirmation of the radiative model. 

Let us now turn to the update of our LO/NLO GRV(94) densities which consists of a 
fine-tuning of the valence-like input densities for f(x, /i 2 ) as well as of the input scale \i. 
This fine-tuning of the input results mainly from new HERA data [7], [| on F 2 (x, Q 2 ) 
and to a certain extent also from new large-x results (such as asymmetry measurements 
of Drell-Yan dilepton production in pp and pal collisions |], [HJ) discussed in more detail 
in Section 3. For the running coupling a s (Q 2 ) at the next-to-leading order (NLO) we 
have utilized the exact numerical solution^] of 

d^m = -^ as{Q) ~mT 2as{Q) (1) 

with /3 = 11 — 2//3 and f3\ = 102 — 38//3, which is nowadays used for NLO analyses 



11, 12, p| since it is stable against yet higher order contributions, thus being more 



appropriate in the low Q 2 region. This is in contrast to using M, 0, H] the approximate 



1 Alternatively, the exact solution can be written implicitly 



i SL ^ 0i i 



47T fa 

0oa s (Q 2 ) 01 



and our new NLO results correspond, for a s (M|) = 0.114, to the following values of A^I : Aj^' ' 6 -* 
299.4, 246, 167.7, 67.8 MeV. 
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solution 

MQ 2 ) 1 _ ft lnhi(Q 2 /A 2 ) 

4vr " fa HQ 2 / A 2 ) (3$ [\n(Q 2 /^ 2 )} 2 1 ] 

which is sufficiently accurate for Q 2 > m 2 . Here, A refers in NLO to A = Aj^g and in 
LO (/3i = 0) to A = Alo We have chosen cn s (M§) = 0.114 for obtaining our exact 
numerical NLO solutions from (1) for Q 2 > ^ LO - This choice, which is slightly preferred 
in our present analysis, agrees with the average value of the space-like momentum-transfer 
measurements [|13|, [L4], [15]], a s (M|) = 0.114 ± 0.005. The statistically dominating time- 
like e + e~ LEP Z°-data imply a somewhat larger 'world average' [fl3| , [TJ[] of cn s (M§) = 
0.118 ± 0.005 with an error which is to some extent uncertain and debatable. It should 
be kept in mind, however, that LEP data (Z° hadronic decays) allow [H| also for a much 
smaller strong coupling, a s (M§) = 0.101 ± 0.013. 

Should significantly higher values of a s , e.g. a s (M|) = 0.118, turn out to be undebated 
and everywhere unique, then our input scale /i will obviously increase closer to 1 GeV 
which compels one to give up the strict valence-like sea input xq(x,fi 2 ), but not the 
valence-like gluon input xg(x,fj, 2 ), in order to reproduce all F%(x, Q 2 ) data at Q 2 > 
1 GeV 2 [|TTJ. Alternatively one may keep the valence-like sea input xq(x,fi 2 ), also for 
a s (M|) = 0.118, as a parameter-free seed for the small-x structure of ^(x, Q 2 > fi 2 ), in 
case one intends to predict and explain the HERA-data merely above Q 2 ~ 3 GeV 2 111 



Furthermore the conventional approximate formula (2), being sufficiently reliable for 

(4,5, 
MS 



Q 2 > m 2 , corresponds in our case to A^' 6 " 1 = 257, 173.4, 68.1 MeV which reproduces the 



exact solutions to even better than 0.5% for Q 2 > 5 GeV 2 . Our LO results correspond to 
A^q 4 ' 5 ' 6 ^ = 204, 175, 132, 66.5 MeV which leads to the (theoretically less relevant) value 
of a^°(Mf ) = 0.125. In both cases we have used for the a s matchings 

m c = 1.4 GeV, m b = 4.5 GeV, m t = 175 GeV. (3) 

These masses are used in all our subsequent LO and NLO analyses for heavy quark 
production. In particular the value of m c , which is slightly lower than the one previously 
employed in [Q, is favored. 

The free parameters of the non-singlet input densities u v , d v , A = d — u and of the 
valence-like singlet input distributions d + u and g at Q 2 = fi 2 have been fixed using the 
following data sets: the published 1994 and 1995 HERA F| results §, § for Q 2 > 2 GeV 2 ; 
the fixed target F 2 P data of SLAC [|TJ, BCDMS [0, |TJ, NMC and E665 |H[ sub- 



ject to the standard cuts Q 2 > 4 GeV 2 and W 2 = Q 2 (- - 1) + m 2 > 10 GeV 2 ; the 
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structure function ratios F£/F p E3, E31, together with the u v /d v results from { vp{d) DIS 



23 1; the Drell-Yan muon-pair production data of E605 |25|] for d 2 a pN j 'dxpdM^+^~ and 



of NA51 pi and E866 110] for the cross section ration a pd /a pp . The input fit parame- 
ters/normalizations of u v and are further constrained by Jq u v dx = 2 and Jq 1 c4,g?x = 1, 
and the ones of the gluon density by the energy-momentum conservation relation 



x u v (x, ji 2 ) + d v (x, /i 2 ) + 2u(x, jj 2 ) + 2d(x, /i 2 ) + g(x, /i 2 ) dx = 1. (4) 

The resulting LO input distributions at Q 2 = /i 2 = 0.26 GeV 2 are then given by 

xu v {x,nl ) = 1.239 x a48 (l -x) 2 ' 72 (1 - 1.8-s/x + 9.5x) 

xd v (x,fx1 ) = 0.614 (1 — x) ' 9 xu v (x, /i 2 ) 

xA(x,/^ ) = 0.23 x a48 (1 - x) 11 ' 3 (1 - 12.0^ + 50.9x) 

x(u + d)(x, / u 2 j0 ) = 1.52 x 015 (l-x) 91 (l - 3.6^ + 7.8x) 

x(?(x,/i 2 ) = 17.47 x L6 (l-x) 3 8 

Ssfo^Lo) = X *( X > Mlo) = ( 5 ) 

where A = J— u. It is interesting to note that /ilo = 2.5 A^b = 2.9 A^b and a^l/^Lo) A" = 
0.24. The corresponding NLO(MS) input at Q 2 = /4 LO = 0.40 GeV 2 is Q 

xu^x, /4lo) = 0.632 x - 43 (1 - x) 3 - 09 (1 + 18.2x) 

XC? t ,(x,/iN LO ) = 0.624 (1 - x) L ° XM t ,(x,/XN LO ) 

xA(x,/4lo) = 0.20 x a43 (l -x) 12 ' 4 (l - 13.3^ + 60.0x) 

x(u + cT)(x,/4lo) = 1.24x a20 (l-a;) 8 ' 5 (l-2.3v / x + 5.7x) 

x<?(x,/4 LO ) = 20.80 X 1 ' 6 (1 - x) 4 - 1 

xs(x,/4 LO ) = xs(x,/i^ LO ) = 0. (6) 

Note again that /^nlo = 2.1 = 2.6 and OL s {lAi^o) = 0-18, and that there is a 

correlation between the chosen value of a s (M|) and the resulting value for //nlo, which 

increases with a g (M§) as already discussed above. Note that our obtained value for /^nlo 

2 The power a soa of the valence-like LO/NLO sea input densities x(w + d)(x, /i 2 ) ~ i""™, as a; — > 0, 
depends strongly on the choice of /i, i.e. on the chosen value for a s (Af§) or, equivalently, on A^ 3 ). 
We would obtain almost equally agreeable results if we continue to use the same a s (Q 2 ) as in our 
previous analyses [[[J |[ [| which was based on the approximate evolution formula (2). Matching the 
approximate NLO(MS) a s (Q 2 ) to the exact numerical solution of eq. (1) at, say, Q 2 = 3 GeV 2 results 
in a s (M§) = 0.110. As discussed above, this a s value lies at the lower end of the presently allowed 
experimental bounds. In this case one obtains /i^LQ — 0-30 GeV 2 and ~ 0.35, instead of 0.20 in 

eq. (6). 
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would have been larger if we had used the (inappropriate) approximate formula (2) in- 
stead of the exact solution of (1). It should be furthermore noted that we have chosen, 
as previously ||, a vanishing strange sea at the input scale /i in order to comply with 
experimental indications [26, 27 of an SU(3)-broken sea. This choice is also supported 



and slightly favored by our input fits and compares well [53] with recent measurements 



2% of s(x, Q >n ). Our s(x, Q ) is thus generated purely dynamically (radiatively) and 
therefore constitutes, for the time being, an absolute, i.e. parameter-free prediction. If 
future experiments may require a finite strange sea input, then our present results for 
s(x, Q 2 ) have to be interpreted as an absolute lower bound for the strange sea. Further- 
more, the charm contribution F 2 (x, Q 2 ) to F 2 is provided by the perturbatively stable 
fixed-order perturbation theory. In LO it derives from the well known photon-gluon 



fusion process p3] 7* g — > cc. For the NLO calculations we employ the O(aj) coefficient 



functions of |30[ as conveniently parametrized in [pl[ . In both cases we use m c = 1.4 GeV 
as given in (3) and choose the factorization and renormalization scales to equal Am 2 . The 
bottom contribution F 2 b is marginal, reaching at most 1 to 2%. 

Our resulting small-x predictions for F 2 are shown in Fig. 3 and for xg (x, Q 2 = 20 
GeV 2 ) by the solid curve in Fig. 2. These results are not too different from our NLO 
GRV(94) expectations when comparing them with the NLO (94) results in Figs. 1 
and 2. Moreover, the NLO results in Fig. 3 are favored over the LO ones in the small 
Q 2 region, Q 2 < 3 GeV 2 (which might indicate that future NNLO contributions could 
even improve our NLO results in the small Q 2 region around 1 GeV 2 ). Furthermore, 
their dependence on choosing different factorization scales [[14], |32[, instead of fi 2 F = Q 2 , is 
obviously weaker than in LO. Nevertheless, the present LO/NLO stability is even better 
than for the GRV(94) F 2 -predictions. The LO and NLO input densities in eqs. (5) and 
(6) are shown in Fig. 4 and compared with the NLO ones of GRV(94) || as well as with 
their evolutions to Q 2 = 5 GeV 2 . The appropriate valence-like NLO input F 2 (x, /i^ixO^ 
which eventually vanishes as x — > 0, is also shown in the lower right corner of Fig. 3. 
This illustrates very clearly that the predictions at x < 10~ 2 and Q 2 > fi^^o are of a 
purely dynamical origin, in particular the increase of F2 with x as x — > is due to the 
non-vanishing input at x > 10~ 2 . Nevertheless, as evident from Fig. 3, our predictions for 
Q 2 < 1 GeV 2 fall below the data in the (very) small-x region. This is not unexpected for 
leading twist-2 results, since nonperturbative (higher twist) contributions to F 2 (x,Q 2 ) 
have eventually to become dominant for decreasing values of Q 2 . 
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It is also interesting to note that the total momentum fractions carried by the NLO 
input distributions at Q 2 = /i 2 in Fig. 4 amount to 56, 30 and 14% for valence, gluon 
and sea densities, respectively, which are very similar in LO and similar to our GRV(94) 
results ||. 

Further typical small-x predictions for xg(x,Q 2 ) and xu(x,Q 2 ) are shown in Fig. 5 
together with their respective inputs at Q 2 = [i 2 which become, particularly for the gluon 
input, vanishingly small at x < 10~ 2 . This illustrates again the purely dynamical origin 
of the small-x structure of gluon and sea quark densities at Q 2 > /i 2 . Also noteworthy 
is the stability of u(x, Q 2 ) at Q 2 ^> /i 2 , i.e. not only the perturbative LO/NLO one 
but also that with respect to our GRV(94) results. This stability is almost as good as 
the one required for a physical quantity like F 2 (x, Q 2 ). The situation is, as usual ||, 
different for g(x, Q 2 ) which is, however, not as relevant since the gluon density is not 
directly measurable. In fact, despite the sizeable difference of the LO and NLO gluon 
distributions in Fig. 5 in the small-x region, the directly measurable gluon-dominated 
heavy quark contribution F%{x, Q 2 ) shows a remarkable perturbative stability even for 
very large values of the (factorization) scale, such as [i 2 F ~ Q 2 , as will be shown in the 
next Section. Furthermore, xg and xu at Q 2 S> /i 2 increase almost linearly for 10~ 5 < x < 
1CT 3 and 1CT 5 < x < 10~ 2 , respectively, on the double-logarithmic plots in Fig. 5. They 
can thus effectively be represented by xf(x,Q 2 ) ~ ^--VO^ 2 ) with effective slopes 

X g (x, 5 GeV 2 ) ~ 0.24 (0.34), X g (x, 20 GeV 2 ) ~ 0.30 (0.39) 

\ u {x, 5 GeV 2 ) ~ 0.21 (0.19), \ a (x, 20 GeV 2 ) ~ 0.26 (0.27) , (7) 

in NLO(LO), valid in the above mentioned x-intervals. These steep (A/ > 0) dynamical 
small-a: predictions are somewhat smaller than the ones of GRV(94) 

Fig. 6 is an alternative way to present our dynamical (steep) small-x predictions using 
the slope of F$(x, Q 2 )- This directly measurable gluon dominated quantity, dF^jd InQ 2 , 
exhibits again a good perturbative stability for Q 2 > 1 GeV 2 . The similarity of our 
new (modified) results with our previous NLO GRV(94) ones is interesting as well. We 
refrain here from plotting any experimental results which are usually Q 2 -averaged |33|, [54]] 



and thus may easily give rise to erroneous and misleading interpretations and conclusions 
when compared with theoretical results for dF%(x, Q 2 )/d InQ 2 which depend strongly on 
the specific choices of x and Q 2 . Nevertheless, our slopes in Fig. 6 are consistent with 
present HERA measurements PBl B41 for not too small values of Q 2 . 
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3 The role of further deep inelastic and hard scattering 
data and very small— x predictions 



As mentioned in the Introduction and the previous Section, fixed-target DIS measure- 
ments and data on F 2 c (:r, Q 2 ), pp(n) — > fi + fi~X and the pp —>■ W ± X charge asymmetry 
are relevant for fixing and testing g(x, Q 2 ), d — u and d/u. In Figs. 7 and 8 we compare 
our fitted LO and NLO results with the fixed-target F$ data[] [17-21] where we show 
only the kinematic region which has been used to determine our valence-like gluon and 
sea (u + d) and the valence input distributions at Q 2 = jj^ Q NLO as given in eqs. (5) and 
(6). The quality of our NLO fits (which are slightly favored over their LO counterparts) 
in Fig. 7 clearly demonstrates that our chosen value of a s (M|) = 0.114 in Section 2 is 
fully consistent with all present fixed-target high precision DIS (non-singlet) data in the 
large-x region. In our present analysis we employ the F^/F^ data |22|, ^] as extracted 
without any (still uncertain) corrections for a possible EMC effect and nuclear shadow- 



ing in deuterium. Together with the d v /u v constraints from pp and vd data this 
procedure leads to a 'traditional' large-x behavior, d v /u v — > as x — > 1. Note that we 
investigated more complicated parametrizations for d v (x,fi 2 ) than used in (5) and (6), 
but found them unnecessary even at the present level of accuracy. 

In Fig. 9 we compare our new (modified) LO and NLO parton distributions with the 
relevant data on DIS charm production []35|, |37|, ^j. From HERA only the most recent 
(preliminary) 1997 ZEUS results |39|] are displayed, since they supersede the previous 
(published) measurements [[37], |38[ by their greatly improved accuracy. As discussed in 
the previous Section, the gluon g(x, ix 2 F ) dominated F^ (x, Q 2 ) is calculated using the fully 
predictive fixed order (LO/NLO) perturbation theory || ^ which also underlies the 
actual analysis of the partial data |37|, ^] utilized for extracting the total F£(x, Q 2 ). The 
factorization scale fi F (being as usual assumed equal to the renormalization scale) should 
be preferably chosen to be \j? F = Am 2 [^J. The resulting predictions in Fig. 9 are in 
perfect agreement with all available data (including the original fixed-target EMC data 
||36|| ) and are furthermore perturbatively stable. Even choosing a very large scale like 
y? F = 4(Q 2 + 4 m 2 ), the NLO results remain essentially unchanged at small-x ||D| as 



3 The normalizations of the F% data sets are allowed to float within their experimental uncertainties; 
the resulting normalization factors are 1.00 (SLAC), 0.98 (BCDMS), 1.01 (NMC, E665). The BCDMS 
data have been taken as analyzed in pg[ , i.e., with a shift of the central values due to the main systematic 
error and correspondingly reduced full errors. This is still the consistent treatment, also in conjunction 
with the new NMC data |l9| . Furthermore we take over the target-mass corrections of |kJ for the SLAC 
and BCDMS Ff data. 
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shown by the 'high scale' dotted curves in Fig. 9. This latter stability renders attempts to 
resum supposedly large logarithms (InQ 2 /m 2 c ) in heavy quark production cross sections 



superfluous [11], [T2|, |41|. It should be noted that the charm (and bottom) production 
data strongly constrain the gluon distribution and will eventually be used to determine 
g(x,fi 2 7 ) directly from experiment ]§[]. 

The asymmetry measurements of Drell-Yan dilepton production in pp and pd collisions 
H |10| which have been instrumental in fixing d — u (or d/u), in particular the very 



recent Fermilab-E866 data |nj for the extented x range 0.03 < x < 0.35, are compared 



with our updated LO and NLO results in Fig. 10. For comparison the consequences of 
our NLO GRV(94) d/u ratio are shown as well which has been originally constrained 
just by the CERN-NA51 measurement || at x = 0.18. Our present new results for 
d/u are also consistent with (although slightly lower than) the recent preliminary semi- 
inclusive HERMES measurements at 0.05 < x < 0.2. The sensitivity of the Drell-Yan 
asymmetry ond-ii can be most easily seen from the LO expression 

_ o pp - o pn 2a pp (u-d)(u-d) + § (uu - dd) /oN 
Ai — — — ' — ' v 1 



a PP _|_ a pn aP d (w + cf) (m + gO + | {uu -dd) + f 

due to a pN oc Y,u,d,s e g [q{ x i)q( x 2) + Q{ x 2)q(xi))- The relevant NLO differential Drell-Yan 
cross section a pN = d 2 a pN /dM^+^-dxp can be found in the Appendix of [43|, except 



for eq. (8a) which has to be modified |14], in order to conform with the usual MS 
convention for the number of gluon polarization states 2(1 — e) in 4 — 2e dimensions. 

Having fixed d — u and d + u as well as the valence densities u v and d v (cf. Figs. 7 
and 8), our strongly constrained u and d distributions are now confronted in Fig. 11 with 
the W ± — > l^v charge asymmetry measurements at the Fermilab pp collider The 
W ± rapidity asymmetry 

= da(e + )/dy e - da(e-)/dy e , . 

[Vi) da(£+)/dy e + da{£-)/dy £ U 

of the charged leptons from the — > l^v decays with the lepton rapidity y^ is a sensitive 
probe of the difference between u and d quark distributions at Q 2 = M^- Our LO and 
NLO predictions are in perfect agreement with present data. A good agreement is also 
obtained [^] with our GRV(94) densities 0], as shown in Fig. 11, and even with our 
previous [u = d) dynamical LO and NLO distributions |], ||. 

Further constraints on the gluon distribution, besides those from DIS and Drell-Yan 



data |P2]|, could be obtained via pp — > jX, pp — > jet + X, etc., but the predicted cross 
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sections are quite sensitive to assumptions concerning the magnitude of the intrinsic 
transverse momentum kx of the partons which is not well understood at present. This 
holds true in particular for the prompt photon production data where, despite the large 
scale (fiR t F) uncertainty [^8], f5H , an additional sizeable /c T -smearing seems to be required 



PD| , |50| , |TT| ] in particular in view of the recent Fermilab (E706) pBe — > 7X measurements 
51 1 . A further interesting source of information, sensitive to g(x, Q 2 ), will be provided by 



the longitudinal structure function Fl(x, Q 2 ) = F2 — 2xF\ as illustrated for example in ||. 
The available data on F^ are unfortunately still of rather limited accuracy || |5^| . Finally 
it should be mentioned that the gluon density is also tested in a NLO determination of 
the strange sea quark density s(x, Q 2 ) from neutrino induced charm production data, 
i.e. opposite-sign dimuon events originating from uN — > fi~cX with c — > /i + z^s. Since 
dimuon events give, among other things, direct access to s(x, Q 2 ) via W + s — > c, W + s — > eg 
and W + g — > cs, etc., they also probe our purely dynamically generated s(x, Q 2 ), i.e. the 
assumed vanishing strange sea input in eqs. (5) and (6), throughout the entire x-region. 
Since our slightly modified input densities in (5) and (6) give rise to s(x, Q 2 ) and g(x, Q 2 ) 
which are similar to the GRV(94) ones, they result in similarly agreeable predictions for 
ViV — > n~ji + X ||28||. It should, however, be kept in mind that a vanishing strange sea 



input is by no means a crucial ingredient of the dynamical approach. 

As our parameter-free small-x predictions for parton distributions at x < 10~ 2 are 
entirely of QCD-dynamical origin and depend, apart from intrinsic theoretical uncertain- 
ties, rather little on the detailed input parameters at x > 10~ 2 , it is interesting to study 
these predictions in kinematic regions not accessible by present DIS experiments and to 
compare them with our previous GRV(94) densities. Of particular interest here is the 
comparison in Fig. 12 at Q 2 = 10 4 GeV 2 and extremely small x, i.e. 10 -8 < x < 10~ 5 , rele- 
vant to questions concerning neutrino astronomy ||53|| . These results at x < 10~ 5 indicate 



that, for example, ultra-high-energy neutrino nucleon cross sections, which are sensitive 
to parton densities at x values as small as 10~ 8 , can be rather reliably calculated to within 
about 20%. This follows not only from the perturbative stability at extremely small val- 
ues of x, where moreover our improved predictions are comparable to the ones based on 
the NLO GRV(94) densities, but also from the fact that the predictions at Q 2 ~ are 
rather independent of the specific choice for the renormalization scale hr appearing in 
a s {n 2 R ) and for the factorization scale \ip appearing in the parton densities f(x,/i F ). We 



have checked this by taking pi R — pL F with [O, |32| Q/2 < /xp < 2Q which requires, of 
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course, also corresponding modifications^ in our input scale \i and /(x,/x 2 ) in eq. (6). 



4 Discussion and Summary 

As demonstrated above, the radiative (dynamical) GRV(94) parton distributions || dis- 
agree with recent precision HERA data only within the margins resulting from the 10% 
uncertainty in their input scale \l 0, [|. Taking into account also the new large-x parton 
constraints and a s -results, we have generated new sets of LO and NLO dynamical parton 
densities corresponding to a s (M z ) = 0.114 for the purpose of future precision analyses. 
Nevertheless one can in practice still utilize the former GRV(94) distributions ||, in 
particular in view of the fact that in most applications of these parton densities to, 
say, high-pr jet, photon or heavy quark production, the usually considered (see, for 
example, [14j |32| for a recent review and comparative discussion) freedom in the choice of 



the factorization and/or renormalization scale, e.g. < hf,r < %Pt, overshadows the 
present modifications of our previous distributions ||. 

It is also necessary to mention that the task of searching the ultimately correct parton 
distributions is not only affected by the above mentioned higher order uncertainties, but 



also by the discrepancies between the data sets used, e.g. between the NMC data |2(J on 



the one hand and the CCFR |55] and HERA @, 0, §] data on the other, which are not 



well understood at present P3L Furthermore, recent attempts to calculate the (non- 
perturbative) input parton densities from first principles using the chiral soliton approach 
yielded, besides the valence densities, also a valence-like sea density in the small-x region 
at Qq = 0.3 — 0.4 GeV 2 - a scale set by the inverse average instanton size [ 57| . It remains 



to be seen, however, whether a sizeable valence-like gluon density at the same 'dynamical' 
input scale is also within the realm of this approach. It is also interesting to remark that 
a valence-like gluon input density, with a momentum fraction compatible with our results 
|T|, |2], |3| , has been obtained from considerations of intrinsic nucleon Fock-states . 

Finally it should be emphasized that the stable parameter-free dynamical predictions 
for parton distributions in the extremely small-x region, 10~ 8 < x < 10~ 5 , allow for 
rather reliable estimates of ultra-high-energy neutrino nucleon cross sections relevant to 
questions concerning neutrino astronomy [53|. 



4 The implementation of this modification amounts to replacing everywhere in Q the NLO a s (Q 2 ) 
and the f(x,Q 2 ) by a s ((Ap) and f(x,n 2 F ), respectively, while the common MS Wilson coefficients in Q 
have to be replaced by C; — * d + InQ 2 / ' \j? F for i — q, g. 
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A FORTRAN package containing our new LO and NLO(MS) parton densities as 
well as F2 ,b (x,Q 2 ), calculated in fixed-order perturbation theory, can be obtained by 
electronic mail on request. Instead of using the appropriate massive quark subprocesses 
for calculating heavy quark production rates in fixed-order perturbation theory, rough 
estimates (valid to within a factor of 2, say) of 'heavy' quark effects can be easier obtained 
with the help of the massless 'heavy' quark distributions c(x, Q 2 ) and b(x, Q 2 ) given in 
For further convenience we also provide the NLO(DIS) distributions which are related to 
the NLO(MS) ones according to eq. (21) of ref. §[. 
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Figure Captions 



Fig. 1 Comparison of our NLO GRV(94) small-x predictions |f§ for the proton structure 
function F 2 with recent precision measurements at HERA || [/]]. For illustration 
we also include the most recent preliminary data [g] (open symbols). The typical 
uncertainties of the GRV(94) predictions are illustrated by the dotted curves which 
are obtained by increasing /znlo by 10%, keeping everything else (valence-like input 
densities, a s , etc.) unchanged. The results of our present new analysis, where all 
these input quantities are consistently modified, lie in between the curves shown at 
small-x. 

Fig. 2 As in Fig. 1, but for the NLO gluon density at Q 2 = 20 GeV 2 . The result of our 
present new analysis is shown as well (solid curve). The shaded bands represent the 
preliminary experimental small-x constraints as extracted from F 2 -measurements 
at HERA and the four data points are derived from deep-inelastic inclusive charm- 
production as analyzed by HI |J. 

Fig. 3 Comparison of our new LO and NLO small-x results for F 2 (x, Q 2 ), arising from the 
inputs (5) and (6), with HERA data for Q 2 > 1 GeV 2 §, g. The valence-like NLO 
input, according to Fig. 4, is shown by the curve (/^nlo) a * the l° wer right corner. 
To ease the graphical representation we have plotted F 2 (x, Q 2 ) +i(Q 2 ) x 0.5, with i 
indicated in the figure. The published data || [7| (closed symbols) are, different from 
Fig. 1, shown with small normalization changes as obtained in our fits: Hl(94)*0.99, 
ZEUS(94)*1.01, shifted vertex data [ZEUS(94) and Hl(95)]*0.97. 

Fig. 4 The valence-like LO and NLO input densities xf (/ = u v , d v , u, d, g) according to 
eqs. (5) and (6) at Q 2 = fx 2 LO = 0.26 GeV 2 and Q 2 = /4 L0 = 0.40 GeV 2 . The strange 
sea s = s vanishes at the input scales Q 2 = A'-lonlo - 

The NLO GRV(94) input || 
is also shown for comparison, as well as the evo luted results at Q 2 = 5 GeV 2 . 

Fig. 5 The small-x behavior of our radiatively generated gluon and sea-quark distributions 
in LO and NLO. The valence-like inputs, according to eqs. (5) and (6) as presented in 
Fig. 4, are shown by the lowest curves referring to /i 2 for illustration. For comparison 
we also show the NLO GRV(94) predictions 0. The results are multiplied by the 
numbers indicated in brackets. 

Fig. 6 The predicted slope dF 2 (x, Q 2 )/d InQ 2 in LO and NLO. The difference between the 
dotted and solid curves is due to the NLO charm contribution dF^(x, Q 2 )/d InQ 2 . 
For comparison we also show the NLO results of GRV(94) 0. The upper part of the 
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figure refers to the values of Q 2 appropriate to the HERA measurements |33|. 53, 11 
and the lower one to a representative fixed value of Q 2 . 



7 Comparison of our LO and NLO fits with fixed-target F p data JT7|, [Tj| F9, [21], [H 
in the Q 2 > 4 GeV 2 , W 2 > 10 GeV 2 region which were used for determining the 
valence and valence-like gluon and sea input in (5) and (6). The data sets are 
shown with their normalization factors as obtained in the fit. Some points with 
large errors, e.g. E665 at large-x, are omitted. 

8 Same as in Fig. 7 but for the measured F^/F^ and d v /u v ratios |22| 24 1. Only 



the NMC results for F^/F^ are shown, as they are much more accurate than the 
corresponding BCDMS and E665 data [ 23fl . The preliminary semi-inclusive HER- 



MES data for d v /u v are displayed as well f35|. For comparison our previous NLO 
GRV(94) results |§ are also shown. 

9 LO and NLO predictions for F£(x, Q 2 ) in fixed-order perturbation theory (7*-gluon 
fusion, etc.) based on our new LO and NLO parton densities using fi 2 F = 4m 2 , 
m c = 1.4 GeV, compared with data from EMC || and ZEUS ||. The NLO 



results based on the NLO GRV(94) parton densities || are close to the present 
LO(NLO) curves at large (small) values of x. The NLO (high scale) curves refer to 
a significantly larger factorization scale, /i^ = 4 (Q 2 + 4m 2 ). 

10 LO and NLO QCD results for a pd /2a pp for A DY in eq. (8) compared with the 



NA51 || and the most recent E866 Jl0| Drell-Yan dimuon production data. The 



experimental (E866) normalization uncertainty of 0.01 has been used in the fit. The 
NLO(94) curves correspond to the d/u GRV(94) ratio || which has been originally 
extracted just from the the NA51 data point at x = 0.18 shown in the right figure. 

11 LO and NLO predictions for the W ± charge asymmetry A(yi) in eq. (9). Note that 
these data have not been used for fixing our new input distributions. The NLO (94) 
predictions, based on the GRV(94) densities ||, are shown for comparison. The 
relevant LO/NLO expressions for A(yn) can be found in ||7j (and references therein). 



The Fermilab-CDF data are taken from [46]. 



12 Predictions for F%(x, Q 2 ) for extremely small values of x. The difference between 
the dotted and solid curves is due to the NLO heavy quark (charm, bottom) con- 
tributions, which derive from photon-gluon (quark) fusion processes. The NLO (94) 
results correspond to the parton densities of 0. The NLO results at these large 
values of Q 2 are insensitive to the specific choice of the factorization scale \ip. 
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